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The methods which are actively used for electronic structure calculations of low- 
lying states of heavy- and superheavy-element compounds are briefly described. The 
advantages and disadvantages of calculations with the Dirac-Coulomb-Breit Hamil- 
tonian, Huzinaga-type pseudopotential, shape-consistent Relativistic Effective Core 
Potential (RECP) and Generalized RECP are discussed. The nonvariational tech- 
nique of the electronic structure restoration in atomic cores after the RECP calcu- 
lation of a molecule is presented. The features of some approaches accounting for 
electron correlation, the configuration interaction and coupled cluster methods, are 
also described. The results of calculations on El 13, El 14, U and other heavy-atom 
systems are presented. 



1. INTRODUCTION 

High-precision calculations of molecules with heavy and superheavy atoms that provide 
"chemical accuracy" (1 kcal/mol or 350 cm~^) for excitation and dissociation energies of 
low-lying states are extremely time-consuming. Employing the latest theoretical and 
program developments is necessary on the following stages: 

(A) selection of an effective spin-dependent Hamiltonian; 

(B) basis set optimization; 

(C) appropriate way of accounting for correlation. 
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In order to minimize the computational efforts necessary to provide a given accuracy 
in calculation of properties, it is important to achieve the equivalent (balanced) level of 
accuracy in each of these stages in the most economical way. Moreover, too high accuracy 
which can be formally attained at the first two stages by, e.g., (a) employing an effective 
Hamiltonian, in which inactive core electrons are treated explicitly or/and (b) using a too 
large basis set etc. can result in abnormal requirements to computers at the last stage. 

In the present paper, the main attention is paid on items (A) and (C). The Dirac- 
Coulomb-Breit (DCB) Hamiltonian and the Relativistic Effective Core Potential (RECP) 
method which are widely employed [at stage (A)] are described in sections 2 and 3. 
The Configuration Interaction (CI) and Coupled Cluster (CC) methods which are most 
popular in correlation calculations, [at stage (C)] are presented in sections 4 and 5. In 
opposite to the density functional approaches, the CI and CC methods allow one to study 
excited electronic states of a given symmetry with high level of accuracy. 

2. DIRAC-COULOMB(-BREIT) HAMILTONIAN 

It is well known that the Dirac-Coulomb (DC) Hamiltonian with the Breit interaction 
and other Quantum ElectroDynamic (QED) corrections taken into account can theo- 
retically provide a very high accuracy of calculations of heavy atoms and heavy-atom 
molecules. The DC Hamiltonian has the form (in atomic units e — m — h — 1, where e 
and m are the electron charge and mass, h is Planck constant): 

H">= = ^h°w+5;i-, (1) 

p p><i 

where indices p, q run over all the electrons in an atom or molecule, rpg is the distance 
between electrons p and g, and the one-electron Dirac operator h° is 

h° = c(a • p) + mc^p - 1) + V''' , (2) 

c is the speed of light, V^'^"'^ is the nuclear potential which can include the effect of finite 
nuclear size etc., p——i'V is the electron momentum operator, a, (3 are the 4x4 Dirac 
matrices. 

The lowest-order QED correction includes the interelectronic exchange by one trans- 
verse photon in the Coulomb gauge and leads to so-called Dirac-Coulomb-Breit Hamilto- 
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nian, 

n^CB ^ ^BC^J-B,,, (3) 

p>q 

where 

ojpq designates the frequency of the photon exchanged between electrons p and q. A 
low- frequency expansion of the cosines yields the incomplete Breit interaction Bpq{^): 

Bpq{0) = -ocp ■ ocg/rpg + ^[otp-ocg- {ocp ■ rpg){cxg ■ rpg)/rlg] /rpg . (5) 

These terms describe the instantaneous magnetostatic interaction and classical retarda- 
tion of the electric interaction between electrons. The contribution from the first term 
(called Gaunt interaction) to transition energies and hyperfine structure (HFS) constants 
can be observed in atomic Dirac-Hartree-Fock (DHF) calculations [1] (Tables 1 and 2). 

The one-electron basis functions in calculations with the DC(B) Hamiltonian are the 
four-component Dirac spinors. The DC(B)-based calculations have the following disad- 
vantages: 

• too many electrons are treated explicitly in heavy-atom systems and too large basis 
set of Gaussian functions is required for accurate description of the large number of 
radial oscillations which valence spinors have in the case of a heavy atom; 

• the necessity to work with the four-component Dirac spinors leads to serious com- 
plication of calculations as compared to the nonrelativistic case. 



3. RELATIVISTIC EFFECTIVE CORE POTENTIALS 

In calculations on heavy-atom molecules, the DC and DCB Hamiltonians are usually 
replaced by an effective Hamiltonian 

Pv Pv^qv 

written only for valence or "valence-extended" (when some outermost core shells are 
treated explicitly) subspace of electrons denoted by indices and gr^; U^^ is an RECP 
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operator simulating, in particular, interactions of the explicitly treated electrons with 
those which are excluded from the RECP calculation. In Eq. (6), 

is the one-electron operator of the nonrelativistic Schrodinger Hamiltonian. Contrary to 
the four-component wave function used in DC(B) calculations, the pseudo-wave function 
in the RECP case can be both two- and one-component. 

3.1. Huzinaga-type potential 

When forming chemical bonds in heavy-atom molecules, states of core electrons are 
practically unchanged. To reduce computational efforts in expensive molecular calcula- 
tions, the "frozen core" approximation is often employed. 

In order to "freeze" core (c) spinors, the energy level shift technique can be applied. 
Following Huzinaga, et al. [2], one should add the effective core operator U^^^ containing 
the Hartree-Fock (HF) field operators, the Coulomb (J) and spin-dependent exchange 
(K) terms, over these core spinors together with the level shift terms to the one-electron 
part of the Hamiltonian: 

Uhuz = (J-K)bncii] + Yl ^"c'i \'PnJj){'Pnclj\ (i-G- Enjj ^ Enjj+Bnjj) , (8) 

where ric, I and j are the principal, orbital momentum and total momentum quantum 
numbers, the Bnjj parameters are at least of order |2£„^/j| and Enjj is the one-electron 
energy of the core spinor (pucij that is frozen. Such nonlocal terms are needed in order 
to prevent collapse of the valence electrons to the frozen core states. As it will be shown 
below, all the terms with the frozen core spinors (the level shift operator and exchange 
interactions) can be transformed to the spin-orbit representation in addition to the spin- 
independent Coulomb term. 

3.2. Shape-consistent radially-local RECPs 

In other RECP versions, the valence spinors are smoothed in the core regions. Consider 
the shape-consistent radially-local (or semi-local) RECP developed by K. Pitzer's group [3, 
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4]. The nodeless numerical pseudospinors ^nvij{f) ^-re constructed of the large components 
fnjj{f) of the valence (y) DHF spinors (one pseudospinor for each I and j): 



fn.ij{r) , r>R, 



^njAr) = < "T'/' . ■ (9) 



'^^E?=o«i'^'' r < Rc : 

where r is the distance between the nucleus and electron. The matching (or core) radius, 
Rc, is chosen near the outermost extremum for the large component and the coefficients 
are taken such that the pseudospinors are normalized, smooth and nodeless. The power 7 
is typically chosen higher than Z + 1 to ensure an efficient ejection of the valence electrons 
from the core region. 

To derive the RECP components Uij, the HF equations are inverted for the valence 
pseudospinors so that (Pnjj become solutions of the nonrelativistic-type HF equations (but 
with j-dependent potentials) for a "pseudoatom" with removed core electrons [5]: 

UiAr) = ^„-i(r) (^^ - ^ + f - J W + K(r) + Snj^j ^n^r) , (10) 

where Z* = Z — Nc, Z is the nuclear charge, Nc is the number of excluded core electrons, 
J and K are the Coulomb and exchange operators on the pseudospinors (fuvij^ ^uyij are 
their one-electron energies (the same as for the original spinors). 

The radially-local RECP operator U^^^ can be written in the form: 

L 1+1/2 j 

VfL^^ + UUr) + Yl E [C//.(r)-C/L.(r)]P,,, P,,= 5] |/jm,>(Zjm,|, (11) 

;=0 j=\l-l/2\ rnj=-j 

where J — L + 1/2, L — l^"'^ + 1 and 1^°^ is the highest orbital momentum of the core 
spinors, mj is the projection of the total momentum. 
Using the identities for the Pij projectors [6]: 



P 



I 



l,j=l±l/2 



21+1 



(/ + ^ ± ^)P^ ± 2Pi I s Pi] , Pi= Yl ■ (12) 

mi=-l 

the RECP operator can be rewritten in the spin-orbit representation, where I and s are 
operators of the orbital and spin momenta, mi is the projection of the orbital momentum. 

Similar to Huzinaga-type potentials, the shape-consistent radially-local RECPs allow 
one to exclude chemically inactive electrons already from the RECP/SCF stage of calcu- 
lations. Moreover, they have the following advantages: 
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1 The oscillations of the explicitly treated spinors are smoothed in the core regions of 
heavy atoms when generating nodeless pseudospinors. Therefore, the number of the 
one-electron Gaussian basis functions may be minimized, thus reducing dramatically 
both the number of two-electron integrals and the computational time. 

2 The small components of the four-component spinors are eliminated and the non- 
relativistic kinetic energy operator is used. The RECP method allows one to use 
a well-developed nonrelativistic technique of calculation whereas relativistic effects 
are taken into account with the help of spin-dependent semi-local potentials. Breit 
and other two-electron QED interactions can be efficiently treated within the one- 
electron RECPs. 

3 In principle, correlations of the explicitly treated electrons with those which are 
excluded from the RECP calculation can be considered within "correlated" RECP 
versions. Reducing the number of explicitly correlated electrons with the help of 
the correlated RECPs is a very promising way to minimize efforts when performing 
high-precision molecular calculations. 

The disadvantages of the semi-local RECPs are: 

1 By now, different versions of the radially-local RECPs provide a comparable level of 
accuracy for the same number of the explicitly treated electrons. It is clear that the 
explicit inclusion of the outer core electrons into the RECP calculation is the way to 
increase the accuracy. However, the extension of the space of these electrons more 
than some limit does not improve the accuracy as is obtained in all our calculations 
with RECPs. The RECP errors still range up to 1000-3000 cm~^ and more even 
for energies of the dissociation of the lowest-lying states and of transition between 
them. 

2 The reliability of the radially-local RECP versions is not high for transitions with 
the excitations in d, /-shells in transition metals, lanthanides, actinides, etc. 

3 Moreover, the direct calculation of such properties as electronic densities near heavy 
nuclei, HFS, and matrix elements of other operators singular on heavy nuclei is 
impossible as a result of smoothing the spinors in the core regions of heavy elements. 
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To overcome the above disadvantages, the Generalized RECP (GRECP) method (see 
subsection 3.3) and the One-Center Restoration (OCR) procedures (see subsection 3.4) 
were developed. 



3.3. Generalized RECP 

It was shown in paper [7] that a requirement for pseudospinors to be nodeless is not nec- 
essary to generate the shape-consistent RECP components. In the case of pseudospinors 
with nodes, the RECP components are singular because division by zero appears in 
Eq. (10). This problem is overcome in the GRECP method by interpolating the potentials 
in the vicinity of these nodes. It was shown both theoretically and computationally that 
the interpolation errors are small enough. This allows one to generate different potentials, 
Unjj and Unjj, for outer core and valence pseudospinors, unlike the conventional RECP 
approach. 

The GRECP operator is written in the form [8]: 

L 1+1/2 
1=0 j=\l-l/2\ 

1+1/2 



+ EE E {[UnJj{r)-UnMr)]PnJj + PnJj[UnJj{r)-UnMr)]} 

Uncij{r) + Un'Jj{r) 



Uc 1=0 j=\l-l/2\ 
L 1+1/2 



~ E E E 

nc,< /=0 j=\l-l/2\ 



- Unjj{r) 



^n'Jj: (13) 



Pncij = E \nJjmj){ncljmj\ . 



The new non-local terms (the second and third lines in the above equation) were added to 
the conventional semi-local RECP operator. These terms take into account the difference 
between the effective potentials acting on the outer core and valence electrons with the 
same I and j quantum numbers. 

The GRECP method allows one to improve accuracy of calculations by regular manner 
when including more outer core shells explicitly into the GRECP calculations. More 
details on the GRECP method can be found in [9, 10]. To compare different effective 
potential versions by accuracy, we carried out both all-electron calculations with the 
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DC Hamiltonian and calculations with RECPs of different groups. The RECP errors in 
reproducing the DHF all-electron results are studied in [9, 10] etc. One can see from our 
atomic HF calculations [11] and correlation calculations on the Hg [12] and Pb [13] atoms, 
that the accuracy of the GRECP is up to an order of magnitude higher than that of the 
other tested RECPs even for the cases when the same number of outermost core shells is 
treated explicitly. 

Results for the eka-thallium atom (E113) are presented in Table 3. The GRECP errors 
are collected into two groups. The errors for transitions without change in the occupation 
number of the 6d shell are rather small. The errors for transitions with change in the 
occupation number of the 6d shell are about 400 cm~^. The latter errors have a systematic 
nature and are connected with the fact that the 60? shell in the present GRECP version is 
described with the help of nodeless pseudospinors. Of course, these errors can be reduced 
significantly if one includes the 5d electrons explicitly in the GRECP calculations. The 
Self-Consistent (SfC) RECP method was suggested in [9, 15], it allows one to minimize 
the above mentioned errors without extension of space of explicitly treated electrons. 
New terms with an operator of the occupation number of the outermost d (or /) shell 
are added to the RECP operator. This method is most optimal for studying compounds 
of transition metals, lanthanides, and actinides. The comparison of accuracy of different 
RECP versions in calculations on the uranium atom can be found in Table 4 and in 
papers [9, 15]. 

A technique for the "Correlated" GRECP (CGRECP) generation was proposed in [9] 
and essential improvements in this technique were made in [10]. The CGRECP for mer- 
cury was generated in the framework of relativistic CC (RCC) method [19]. The GRECP 
components were constructed for 5s, 5p, 5d, 6s, 6p, 6d, 5/, 5g electrons of the Hg atom. 
The 5s, 5p pseudospinors are "frozen" in calculations with this CGRECP (in fact, they 
are completely excluded from the calculations with the help of the level shift technique 
(see section 3.1 and Refs. [9, 20])) and only 12 external electrons of the Hg atom should 
be explicitly correlated instead of 34 ones in the case of the DC calculations because the 
correlations for the 4/, 5s, 5p electrons and between the 4/, 5s, 5p and 5d, 6s, 6p electrons 
are taken into account by the "correlated" GRECP. It allows one to reduce drastically 
the computational efforts necessary for the "chemical" accuracy in calculations of mer- 
cury compounds. Results of our test calculations are presented in Table 5. One can 
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see that energies of the one-electron excitations in calculations of the Hg atom with the 
all-electron DC Hamitonian for 34 external electrons correlated explicitly by the RCC 
method (DC/34e-RCC) can be reproduced with the accuracy withing 270 cm~^ in the 
12e-RCC calculations when using the present CGRECP. Contribution of the core corre- 
lation effects can be seen from comparison of the DC/12e-RCC and DC/34e-RCC results 
or of the GRECP and CGRECP results. 

3.4. Nonvariational One-Center Restoration of electronic structure in cores of 

heavy-atoms in a molecule (NOCR) 

In the valence region, the electronic density obtained from the two-component GRECP 
(pseudo)wave function very accurately reproduces the corresponding all-electron four- 
component density. In the core region, the pseudospinors are smoothed, so that the 
electronic density with the (pseudo)wave function is not correct. 

The following restoration scheme was developed (see [21, 22] and references): 

• Generation of equivalent basis sets of atomic (one-center) four-component spinors 



parts of the atomic Dirac spinors and l'—2j—l) and two-component pseudospinors 



component GRECP/HF calculations of the same valence configurations of the atom 
and its ions. 

• The molecular pseudospinorbitals (pi are then expanded in the basis set of the one- 
center two-component atomic pseudospinors (for r<R^^^^, where R^^^^>Rc), 

Lmax ^ + 1/2 

0i(x) J2 Yl 'iljmjfnlj{r)Xljmj , (14) 

1=0 j=\l-l/2\n,mj 

where x denotes spatial and spin variables. 

• Finally, the atomic two-component pseudospinors are replaced by the equivalent 
four-component spinors in the molecular basis and the expansion coefficients c^^ij^. 




(where fnij, Qnij are the radial parts, Xijmj are the spin-angular 



{fnij{'''')xijmj} by atomic finite-difference (numerical) all-electron DHF and two- 
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from Eq. (14) are preserved: 



1=0 j=\l-l/2\n,mj 



1+1/2 

'nljnij 



fnlj{^)Xljmj 
\9nlj{r)Xl'jm^ 



(15) 



The molecular four-component spinors constructed this way are orthogonal to the 
inner core spinors of the heavy atom, as the atomic basis functions used in Eq. (15) are 
generated with the inner core electrons treated as frozen. The properties described by the 
operators singular close to (heavy) nuclei are calculated with the restored bispinors 0^ . 
More advanced technique of the variational restoration is proposed in [22]. 



4. CONFIGURATION INTERACTION 

The many-electron wavefunction ^'^^ in the CI method is presented by a linear com- 
bination of determinants Dj 

^ci^j-Cf'D,, (16) 

are some numbers (CI coefficients). In turn, each determinant is an anti-symmetric 
production of one-electron basis functions where N is the number of electrons in the 
considered system. The CI equations are written as 

J2HijCT = E^'Cf', (17) 
J 

where Hjj are Hamiltonian matrix elements in the basis set of the determinants and E^^ 
is the CI energy. To find the coefficients and the energy in the CI method, one should 
diagonalize the Hamiltonian matrix. 

If all the possible determinants are considered then the method (called Full-CI) will 
provide the "exact" solution in the framework of a given one-electron basis set and an 
employed Hamiltonian. However, requirements to the computational resources in the 
Full-CI case are usually so huge that such calculations are practically impossible for 
systems of interest except the cases of very small numbers of correlated electrons and 
basis functions. In almost all the CI calculations, only some selected (the most important) 
determinants are explicitly considered. To take into account the effect of the unselected 
determinants, various semi-empirical corrections (e.g., the Davidson correction [23]) can 
be employed. In precise calculations, the number of selected determinants reaches a few 
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millions and more, therefore a very large Hamiltonian matrix should be diagonalised. The 
iterative diagonalization (Davidson) method is then used to obtain a few low-lying roots 
of this matrix. 

There are two main categories of the CI method [24]: 

• "Conventional CI": the Hamiltonian matrix elements are calculated once and saved 
in memory, 

• "Direct CI": only those Hamiltonian matrix elements are calculated at each step of 
the diagonalization procedure which are required at the moment. 

The CI method has the following advantages: 

1 simplicity of the method, solutions are always exist independently of the number of 
open shells; 

2 it well describes "static" (avoided crossing of terms) and "nondynamic" electron 
correlations. 

The disadvantages of the CI method are: 

1 it is badly working for large number of correlated electrons (when semi-empirical 
corrections on unselected determinants are large); 

2 unsmoothness of potential curves is a result of selection of determinants by some 
thresholds; 

3 the above semi-empirical energy corrections cannot be used when calculating other 
than spectroscopic properties. 

5. THE COUPLED-CLUSTER APPROACHES 

The complete space of {Dj} is divided into two subspaces: 

M.0 , model space, consists of small number (M) of the most important determinants 
{Dm}m=i to describe static and nondynamic correlations, which are taken into 
account exactly on Aio; 
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A4q , rest of space (usually very large), is included approximately to account for dynamic 
correlations (i.e. correlations at small interelectronic distances, "Coulomb holes"). 

The eigenstates of interest are presented as 

M 

l^cc^ =^C„exp[r(-)]|i^^) , (18) 

m=l 

where T^"*) = Ti'^^+T^-^^+ ... is the cluster operator: 

(19) 

ij,ab 

where a+ and cij are the creation and annihilation operators (their combination a^cij 
will replace the i-th one-electron state in the determinant by the a-th one). The coeffi- 
cients {tl^\t\^lh}, etc. are called the cluster amplitudes and are calculated solving Bloch 
equations: 

M 

UHU = HU , (U = Vexp[T(™)]|Dj(D„|) . (20) 



m=l 



The coefficients Cm and final energy are obtained from diagonalization of some 
effective Hamiltonian « on the model space: 

M M 

J2 CJDJ = J2 C,r^\Drn) , (Hf„ = (D,|(exp[-T(-)]Hexp[T(™)])|D^)) . 

m=l m=l 

(21) 

If all the Tjf^'^ are considered in the T^"*) operator then the CC method is equivalent 
to the Full-CI one. However, in practical calculations, the third and following terms in 
j'irn) (three-body and higher order cluster amplitudes) are usually neglected. Such a CC 
version is called CC-SD. There are three basic CC categories [25]: 

• One-state or state-selective; 

• Fock-space or valence universal methods; 

• Hilbert-space or state-universal approaches. 
The CC method has the following advantages: 
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1 It is the size-extensive method, i.e. the energy of the system is scaled properly with 
increase in the number of electrons (whereas the CI method is not size-extensive in 
a general case). 

2 The CC-SD method takes into account the contributions not only from the deter- 
minants obtained from the model space by applying the (l-l-T^^^^^-l-Tg"^^) operator 
but also approximately from all the rest determinants (whereas the CI method with 
the same number of unknown coefficients does not) . 

3 The CC method is one of the best methods for accounting the dynamic correlation. 
The disadvantages of the CC method are: 

1 This is a nonvariational method, i.e. the CC energy is not an upper bound to the 
exact energy of the system (whereas the CI energy is). 

2 The CC equations are nonlinear and the effective Hamiltonian is non-Hermitian. 

3 Intruder states (i.e. such states from the A4q subspace, which are lying within the 

A4o subspace energy span) destroy the convergence of the CC iterations. Alleviation 

the problem is in using: 

— Incomplete model space procedures; 

- Energy shifting, RLE [26], DIIS [27, 28], IPM [29] procedures. 

6. SOME PRACTICAL CALCULATIONS 

Calculations of the spectroscopic constants for the ground and lowest excited states of 
the HgH molecule and for the ground state of the HgH+ ion were carried out with the help 
of the GRECP and RCC [19] methods in [30]. The results are within a few mbohr from 
the experimental data for bond lengths, tens of wave numbers for excitation energies and 
vibrational frequencies. It is demonstrated that the triple cluster amplitudes for the 13 
outermost electrons and corrections for the Basis Set Superposition Errors (ESSE) [31, 32] 
are necessary to obtain accurate results for this molecule. The accurate GRECP/CI 
calculations of the spectroscopic constants for the ground state of the TIH molecule are 
presented in [20], in which the reliability of the semi-empirical energy corrections is in 
particular investigated. 
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The NOCR scheme was apphed in the GRECP/RCC calculations of the P,T-odd 
properties for the TIF molecule [21]. The corresponding GRECP/HF/NOCR results 
are in good agreement with the all-electron DHF results of other groups. Inclusion of 
electron correlation has changed the values on 20%. The previous XOCR version was 
employed in the GRECP calculations of the P, T-odd parameters and HFS constants for 
the YbF [33, 34] and BaF [35] molecules. A reasonable agreement with the experimental 
data for the HFS constants was attained. It was demonstrated that the spin-correlation 
effects of the unpaired electron with the deeply-lying outer core 5s and 5p shells of heavy 
atom should be taken into account in order to perform accurate calculations of the HFS 
and P, T-odd constants. 
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Table 1. Transition energies of the Tin {Z—50), Lead {Z=S2) and Eka-lead (Z=114) atoms 
calculated by the DHF method with Coulomb and Coulomb-Gaunt two-electrons interaction for states 
with the ns^np^ configuration (in cm~^). 



configuration J DC DCG absolute relative (%) 

difference difference 
Tin 





3153 


40 


1 3 

J- .«J 


(5s?/25/>}/25P3/2) 1 











i^4/2^Pl/2^Pl/2) 2 5143 


5139 


-4 


-0.1 


(5s?/25P3/2) 2 5941 


5893 


-48 


-0.8 


(5^V25p1/2) 15873 


15820 


-53 


-0.3 




Lead 




iGsl^M/2) 











(6.2/26p}/26p^/2) 1 4752 


4644 


-108 


-2.3 


(6.2^26^1/26^3/2)2 9625 


9514 


-111 


-1.2 


(6.2^26^3/2) 2 18826 


18592 


-234 


-1.2 


(6.2^26^3/2) 28239 


27995 


-244 


-0.9 



Eka-lead 



(7.2/27p2/2) 














(7.^/27^1/27^3/2) 


1 27198 


26806 


-392 


-1.4 


(7.2/27pJ/27pl/2) 


2 30775 


30391 


-384 


-1.2 


(7s?/27p1/2) 


2 66068 


65225 


-843 


-1.3 


(7s?/27p1/2) 


74527 


73674 


-853 


-1.1 
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Table 2. HFS constants in the Indium (Z=49), Thallium (Z=81) and Eka-thallium (Z=113) atoms 
calculated by the DHF method with Coulomb and Coulomb-Gaunt interaction for different 



configurations (in MHz). 

configuration DC DCG absolute relative (%) 

difference difference 
Indium 

{5sl/^5p{/^) 1913 1900 ^13 ^0?f 

(5s2/2^/2) 288 287 -1 -0.3 

(5s2/2^/2) 4.41 4.40 -0.01 -0.2 

(^^1/2^/2) 1-88 1-88 0.0 0.0 

(5s2^26sv2) 1013 1011 -2 -0.2 

Thallium 

(6s2^26p}/2) 18918 18691 -227 -1.2 

(65^^26^3/2) 1403 1391 -12 -0.9 

(6^2^264/2) 20.8 20.8 0.0 0.0 

(652^264/2) 8.72 8.70 -0.02 -0.2 

(652^2^*1/2) ''826 7807 -19 -0.2 

Eka-thallium 

(7s2/2'^^'}/2) 150168 147538 -2630 -1.8 

(7s2^2''^'3/2) 2007 1983 -24 -1.2 

(752/2^4/2) 34.3 34.2 -0.1 -0.3 

(7s2/2''4/2) 13-5 13.5 0.0 0.0 

(752,285!) 28580 28473 -107 -0.4 



"The magnetic moment /uat and spin / for the Eka-thallium nucleus were taken as those for Thallium. 
The presented results can be easily recalculated as only the proper values of JjLn and / are known because 
they just include the iin/I coefficient. 
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Table 3. Transition energies between low-lying configurations of the eka-thallium (E113) atom 
derived from all-electron calculations and the errors of their reproducing in calculations with different 
RECP versions. All values are in cm~^. 











21e- 




All-el. 


Gaunt 


21e- 


RECP 




DHFG" 


contrib. 


GRECP'' 


of Nash 

et al." 


Configuration 


Transition 


DHFG 


Absolute errors 




energies 


- DHF 






64/2K/27'^l/27^'l/2('^ = 1/2) ^ 










c j4 a^fi 7„2 7„1 ( j o lrj\ 

0^/20^5/2 '^1/2 %/2i'^ - 




Q /I 7 

o4/ 


-23 


282 


64/2K/27Sl/28Sl/2('^ = 1/2) 


34962 


374 





-186 




50316 


395 


6 


148 


Dctg ^2 oa^ ^2 ' ^ 1 /2 [nonrel .av.) 


52790 


395 


6 


148 


Qd^l^M^I^T s^l^Td} {nonrel.av.^ 


4ozio 


o95 


6 


161 


6d3/2 6^5/2 7s (-^ - ^) 


r7i on 


095 


6 


148 


Rrl^ (3 j6 7„1 7„2 / 7 1 /r)\ 

^"3/2""5/2'^l/2'^'l/2(.'^ - i/2j 


f:i /loo 


-DU 


32 


4830 


6^/2 6^/2 7s 1/2 7^1/2 7^3/2 (^eL aw.) 


QQ 1 77 


/I Q 

/!4o 


-4 


5177 


6(^3/2 60(5/2 7s 1/2 7^3/2 M. aw.) 




d24 


-9 


5729 


64/264/27sJ/27pJ/2 (reLaw.) 


iiO /4U 


ZOO 


-2 


5161 


64/26^5/2751/27^3/2 (reL aw.) 


149526 


678 


-10 


5811 


64/2K/2 7si/2(o' = 1/2) 


2d4do5 


796 


-4 


6151 


cj5 7„2 7„2 (' T K /0~\ 

^"3/2""5/2'^l/2'?'l/2W - ^J/^) 


4 / 41U 


778 


403 


-2389 


0"3/2""5/2 ' ^1/2 '^'1/2 iP3/2[rel.av.) 


7/1 OQO 


A OA 
-4Z4 


341 


-2089 


6«3/2 64/2 7s 1/2 7^3/2 M-«^-) 


1104,35 


-6 


306 


-1556 


n j3 a j6 '7 2 7 2 / 7 /o\ 

64/264/27s^/27p^/2(J = 3/2) 


78862 


-416 


375 


-2272 


fi/-73 ^i,-76 '7ti2 7;>^1 (f^^l rtri, \ 
0"3/2""5/27Sl/27Pl/27P3/2l^e/.aW.) 


1 O/i no7 

iU4uy / 


-OD 


405 


-1968 


64/264/27s2/27p2/2(reLaw.) 


137083 


306 


473 


-1436 


64/264/27s?/27Pv2M-a^-) 


110139 


-407 


380 


-2317 


64/264/27«?/27P3/2M.aw.) 


150116 


45 


338 


-1679 


64/2 64/2 7.s2/27p}/2(reZ.aw.) 


139841 


-65 


439 


-2184 


64/264 /2 7s2/27p^/2 (reZ.aw.) 


177157 


361 


506 


-1541 


64/264/27s2/2(J = 5/2) 


239509 


158 


408 


-1603 


64/264/27si/2(-^ = 3/2) 


267208 


481 


579 


-1431 



"All-clcctron Dirac-Hartree-Fock-Gaunt (DHFG) calculation with Fermi nuclear charge distrtibution 
for A = 297. 

''GRECP generated in the present work from DHFG calculation. 

'^RECP from [14] (generated from DHF calculation without Gaunt iteraction). 
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Table 4. Transition energies between states of U (averaged over nonrelativistic configurations) 
derived from all-electron DHF calculations and the errors of their reproducing in calculations with 
different RECP versions. All values are in cm~^. 







RECP of 


Energy- 




Quadratic 


"Frozen 




DHF 


Ermler 


adjusted 


SfC 


SfC 


core" 






et al. [16] 


ppa 


GRECP 


GRECP 


if) 


if) 


Num. of el-ns 


All 


14 


32 


24 


24 


24 


24 


Conf 

Will . 


Tr pnprcrv 


A 


Dsolute error 


5f37o2g(7l _^ 
















5f^7s'^7p^ 


7383 

J <JCJ<J 


00 / 


-4yo 


-00 


-00 


o 
z 


1 A 

14 


UJ 1 o 


36159 


ooz 


lOU 


A 
4 


R 
D 


Q 



1 R 
ID 




1 39QQ 


-lyz 


1 c;/l 


o 
-o 


-0 


-1 


1 R 
-ID 


Kf37„l«^l7„l 

(J / to \jfJj I jJ 


1 798Q 


1/1/1 


KOI 
-OZl 


Q1 


Q1 
-ol 


-1 


-0 


5 f^QcP 


54892 


-izi 


-oyo 


1 /I 
-14 


1 

-10 


1 


01 
-Zl 


^/■37e2R^l . 

O J to \JUj ^ 


















16483 


1 7fi 


7SS 


79 
- 1 zo 


n 
u 


04 


1 S7 
10 t 


5 /•47s2 _^ 


















15132 


- 1 Oo 


87 


1 1 

1 1 


1 1 

- 1 1 


1 R 
-ID 


-oO 




15016 


Qn 
yu 


-44o 


"^7 


-ZD 


1 

-1 



-z 


(J J \JiJj 


34099 


1 0fi7 

-izo / 


-IDo 


OS 
Zo 


1 Q 
-10 


-ZD 


RO 
-DZ 


5f"^6(f7n^ 


32341 


70/1 

- / y4 


/1 1^7 
-40 / 


1 1 
-11 


-ZO 


1 7 
-1 / 


-oy 


OJ i o xjLL ^ 


















3774 


ouy u 


74S 

- 1 40 


1 7 


1 7 


QO 
yu 


QR 
-y D 


5f27s2g^2 ^ 
















5 f27s2«^i7„i 


12646 


-441 


-626 


-16 


-15 


-5 







49638 


-498 


155 


24 


25 


-5 


1 


Kf27„l«^3 


1 06Q7 


608 


-240 


-10 


-10 


13 


1 




19319 


390 


-826 


-26 


-26 


6 





bf7s^Qd^ 
















5f7s^6d^ 


29597 


11666 


-1526 


-896 


-104 


466 


48 


5/1 75^6(^3 ^ 
















5/i7s26d27pi 


18141 


-1367 


-778 


46 


49 


-2 


-2 


5/1 7^26^2 


49158 


-1355 


173 


70 


73 


-3 


-2 


5/1 7s led'' 


7584 


1655 


-331 


-39 


-40 


22 


14 


5/i7si6(i37pi 


21154 


779 


-1055 


-11 


-11 


16 


10 


5/375^6^1 ^ 
















5/5 


100840 


430 


1453 


-1860 


22 


105 


291 



"PseudoPotential (PP) from [17] (generated from all-electron calculation in the framework of Wood- 
Boring [18] approximation). 
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Table 5. Transition energies from all-electron DC and different GRECP calculations of the 
lowest-lying states of the mercury atom and its ions in the [7, 9, 8, 6, 7, 7] correlation basis set from [12] 
for the 12 and 34 correlated electrons" by the RCC method. All values are in cm^^. 



State (leading 


DC 


DC 


GRECP CGRECP 


conf. , term) 


34e-RCC 12e-RCC 12e-RCC 


12e-RCC 


OCt OS Do J — ^ 
00, OS op [ Jr^o ) 










37471 


37208 


37244 


37742 


5(f^6s^6p^^Pi) 


39318 


38992 


39025 


39573 


5d^°6s^6p^{^P2) 


44209 


43675 


43710 


44453 


5d^^6s^6p^^Pi) 


55419 


54769 


54780 


55466 


5d^Hs\'^Sy2) 


84550 


83885 


83919 


84774 


5d'Hs\^S,/2) ^ 










5di06pi(^A/2) 


52025 


51515 


51559 


52059 


5di06pi('^3/2) 


61269 


60476 


60532 


61320 


5diO(i5o) 


151219 


150132 


150202 


151262 



"This number is smaller by one or two for Hg"*" or Hg^+ ions, respectively. 



